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^ : ABSTRACT 

^ I We use Owens Valley Radio Observatory (OVRO) cosmic microwave background 

' (CMB) anisotropy data to constrain cosmological parameters. We account for the 

' OVRO beamwidth and calibration uncertainties, as well as the uncertainty induced by 

■ the removal of non-CMB foreground contamination. We consider open and spatially- 
T-H ■ flat- A cold dark matter cosmogonies, with nonrelativistic-mass density parameter Q,q in 

, the range 0.1-1, baryonic-mass density parameter fi^ in the range (0. 005-0. 029)/i~^, 

! and age of the universe to in the range (10-20) Gyr. Marginalizing over all parameters 

^ I but f^O) the OVRO data favors an open (spatially-flat-A) model with Jig — 0.33 (0.1). 

' At the 2 a confidence level model normalizations deduced from the OVRO data are 

■ mostly consistent with those deduced from the DMR, UCSB South Pole 1994, Python 

'• I-III, ARGO, MAX 4 and 5, White Dish, and SuZIE data sets. 
^ ■ 

I ' Subject headings: cosmic microwave background — cosmology: observations — large-scale 

■ structure of the universe 

> ■ 

"-Hi 1. Introduction 

■ Cosmic microwave background (CMB) anisotropy measurements have begun to provide inter- 
esting constraints on cosmological parameters.'' Ganga et al. (1997a, hereafter GRGS) developed a 
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technique to account for uncertainties, such as those in the beamwidth and the cahbration, in hkeh- 
hood analyses of CMB anisotropy data. This technique has been used with theoretically-predicted 
CMB anisotropy spectra in analyses of the Gundersen et al. (1995) UCSB South Pole 1994 data, 
the Church et al. (1997) SuZIE data, the Lim et al. (1996) MAX 4+5 data, the Tucker et al. (1993) 
White Dish data, the de Bernardis et al. (1994) ARGO data, and the Piatt et al. (1997) Python 
I-III data (GRGS; Ganga et al. 1997b, 1998; Ratra et al. 1998, 1999a, hereafter R99a; Rocha et 
al. 1999, hereafter R99). A combined analysis of all these data sets, excluding the Python data, is 
presented in Ratra et al. (1999b, hereafter R99b). 

In this paper we present a similar analysis of CMB anisotropy data from the OVRO obser- 
vations (Leitch et al. 2000, hereafter LOO). The OVRO detectors and telescopes are described in 
Leitch (1998) and LOO; here we review information about the experiment that is needed for our 
analysis. 

OVRO data were taken in two frequency bands, one centered at 14.5 GHz (Ku band), the 
other at 31.7 GHz (Ka band). Thirty-six fields, along an approximate circle at declination 6 ~ 
88° centered on the North Celestial Pole (NCP) were observed. In our computations we use the 
coordinates for the 36 fields given in Table 2 of LOO. The OVRO measurements were made by 
switching the beam in a two-point pattern along the circle, resulting in a three-beam response to 
the sky signal. The beamthrow is 22'. 16. The zero-lag window function parameters for the OVRO 
experiment are given in Table 1. This and other window functions are shown in Fig. 18 of LOO. 

LOO use multiepoch VLA observations to detect and remove non-CMB discrete source contam- 
ination from the OVRO data. We have also analyzed the OVRO data ignoring 3 of the 36 fields 
that were affected by the strongest variable discrete source; cosmological constraints derived from 
this restricted OVRO CMB anisotropy data set are very consistent with those derived from the 
full OVRO CMB anisotropy data set, so we do not discuss this restricted OVRO data set analysis 
further. 

Since OVRO data were taken at two frequencies, it is possible to fit the data to both a non- 
CMB foreground component (parametrized by the frequency dependent temperature anisotropy 
ATfore oc u'^) and a CMB anisotropy component with spectral index P = 0.® We use the method 
in § 11 of LOO to extract the CMB anisotropy component in the OVRO data, marginalizing over a 
foreground spectral index in the range —3 < /3 < 2 in our likelihood analysis. ^'^'^ 



*See LOO and Mukherjee et al. (2002) for discussions of foreground contaminants in the OVRO microwave data. 

® Although the data themselves are unable to rule out more negative values of beta (LOO, Fig. 14), Leitch et al. 
(1997) use low frequency maps of the NCP region to rule out such values. 

^"Following Mukherjee ot al. (2002) wc have also analyzed the 31.7 GHz OVRO CMB anisotropy data while 
marginalizing over possible 100 /nm and 12 /nm foreground contaminant template (Schlegel, Finkbeiner, & Davis 
1998) correlated components. The cosmological constraints from these analyzes are quite consistent with results 
presented here. This is because although the foreground signal inferred in our analysis is not entirely fit by the 
dust data, they are significantly correlated, and the 31.7 GHz data, modelled either way, is almost entirely CMB 
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CMB anisotropy constraints are derived from the foreground-corrected 31.7 GHz data. The 
31.7 GHz beam profile is weh approximated by a circular Gaussian of FWHM 7'.37 it 0'.26 (one 
standard deviation uncertainty). We use the method of GRGS to account for the OVRO beam 
uncertainty. 

As discussed in LOO, the noise in the 31.7 GHz data indicates the presence of a component 
that is correlated between neighboring fields (this component is small compared to the uncorrelated 
noise in a single scan of data). As a result the 31.7 GHz OVRO data show only one-half of the 
anticorrelation for nearest neighbor fields expected for a triple beam chopped experiment. This 
one-offdiagonal correlated noise is included and its amplitude marginalized over in our analysis. 

A constant offset is removed from the OVRO data; we marginalize over the amplitude of the 
offset to account for this in our likelihood analysis. The 1 cr absolute calibration uncertainty of the 
OVRO data is 4.3%, and the method developed by GRGS is used to account for it. 

In § 2 we summarize the computational techniques used in our analysis. See GRGS and R99a 
for detailed discussions. Results are presented and discussed in § 3. We conclude in §4. 

2. Summary of Computation 

In this paper wc focus on a spatially-flat CDM model with a cosmological constant A.^^ As a 
foil we also consider a spatially open model with no A (see, e.g., Gott 1982; Ratra &; Peebles 1995). 
These models are discussed in more detail in R99a, R99b, and R99. 

The CMB anisotropy spectra in these models are generated from quantum fluctuations in 
weakly coupled fields during an early epoch of infiation and so are Gausssian (sec, e.g., Ratra 
1985; Fischler, Ratra, & Susskind 1985). Consistent with this, the observed smaller-scale CMB 
anisotropy appears to be Gaussian (sec, e.g., Park et al. 2001, Wu et al. 2001; Shandarin et al. 
2002; Polenta et al. 2002), and the experimental noise also appears to be Gaussian, thus validating 



anisotropy. The OVRO data at its two frequencies are shown in Fig. 13 of LOO, the deduced CMB anisotropy and 
foreground signals axe shown in Fig. 16 of LOO, and the dust-correlated emission is shown in Fig. 1 of Mukherjee et 
al. (2002). 

^^At fi = —2.2 for the foreground contaminant, 96% of the 31.7 GHz data is CMB anisotropy. 

^^Models that neglect these correlations are grossly discrepant with the data, while when these correlations are 
accounted for the model fits are reasonable and consistent with the data. This can be seen from Fig. 19 of LOO and 

we find the same. 

"See, e.g., Peebles (1984), Efstathiou, Sutherland, & Maddox (1990), Stompor, Gorski, & Banday (1995), Ratra 
et al. (1997), Sahni & Starobinsky (2000), Carroll (2001), and Peebles & Ratra (2002). While not considered in this 
paper, a time- variable dark energy dominated spatially-flat model is also largely consistent with current observations 
(see, e.g., Peebles & Ratra 1988; Ratra & Quillen 1992; Steinhardt 1999; Brax, Martin, & Riazuelo 2000; Huterer & 
Turner 2001; Chen & Ratra 2002; Deustua et al. 2002). 
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our use of the GRGS likelihood analysis method. 

As discussed in R99a, the spectra are parameterized by their quadrupole-moment amplitude 
Qrms-PSi the nonrelativistic-mass density parameter Qq, the baryonic-mass density parameter Qb, 
and the age of the universe to. The spectra are computed for a range of spanning the interval 
0.1 to 1 in steps of 0.1, for a range of Q.Bh^ [the Hubble parameter h = i7o/(100 km s^"*^ Mpe^"*^)] 
spanning the interval 0.005 to 0.029 in steps of 0.004, and for a range of to spanning the interval 
10 to 20 Gyr in steps of 2 Gyr. In total 798 spectra were computed to cover the cosmological- 
parameter spaces of the open and fiat-A models. Examples of spectra are shown in Fig. 2 of R99a, 
Fig. 1 of R99b, and Fig. 2 of R99. 

Following GRGS, for each of the 798 spectra considered the "bare" likelihood function is com- 
puted at the nominal beamwidth and calibration, as well as at a number of other values of the 
beamwidth and calibration determined from the measurement uncertainties. The likelihood func- 
tion used in the derivation of the central values and limits is determined by integrating (marginal- 
izing) the bare likelihood function over the beamwidth and calibration uncertainties with weights 
determined by the measured probability distribution functions of the beamwidth and the calibra- 
tion. See GRGS for a more detailed discussion. The likelihoods are a function of four parameters 
mentioned above: Qrms-PS; ^o-, ^Bh^, and tg. We also compute marginalized likelihood functions 
by integrating over one or more of these parameters after assuming a uniform prior in the relevant 
parameters. The prior is set to zero outside the ranges considered for the parameters. GRGS 
and R99a describe the prescription used to determine central values and limits from the likelihood 
functions. In what follows we consider 1, 2, and 3 a highest posterior density limits which include 
68.3, 95.4, and 99.7% of the area. 



3. Results and Discussion 

Table 2 lists the derived values of Qrms-PS and bandtemperature STi for the fiat bandpower 
spectrum, for the OVRO data. These numerical values account for the correlated noise and offset 

removal, the beamwidth and calibration uncertainties, and the uncertainty due to non-CMB diffuse 
foreground contamination removal. These results are very consistent with those of LOO. For the flat 
bandpower spectrum the OVRO data average 1 a 6Ti error bar is ~ 14%^^ : OVRO data results 
in a very significant detection of CMB anisotropy, even after accounting for the uncertainties listed 
above. 

As discussed in R99a, R99b, and R99, the four-dimensional posterior probability density distri- 
bution function -Z^(Qrms-PS) ^o,^Bh'^,to) is nicely peaked in the Qrms-PS direction but fairly flat in 
the other three directions. Marginalizing over Qrms-PS results in a three-dimensional posterior dis- 



^■^For comparison, the corresponding 1 a STi error bar is ~ 10 — 12% for DMR (depending on model, Gorski et al. 
1998), ~ 15% for ARGO (R99a), and - 14% for MAX 4+5 (Ganga et al. 1998). 
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tribution L{i),Q, ^sh"^, to) which is steeper, but still relatively flat. As a consequence, limits derived 

from the four- and three-dimensional posterior distributions are generally not highly statistically 
significant. Wc therefore do not show contour plots of these functions here. Marginalizing over 
Qrms-PS and one other parameter results in two-dimensional posterior probability distributions 
which are more peaked (see Figs. 1). As in the ARGO (R99a), Python (R99), and combina- 
tion (R99b) data set analyses, in some cases these peaks are at an edge of the parameter range 
considered. 

Figure 1 shows that the two-dimensional posterior distributions allow one to distinguish be- 
tween different regions of parameter space at a fairly high formal level of confldence.^^ For instance, 
the open model near Qq ~ 0.75, Qsh'^ ^ 0.03, and to ~ 20 Gyr, and the flat-A model near Qq ~ 0.6, 
O^sh^ 0.03, and to ~ 20 Gyr, are both formally ruled out at 3 cj confidence. However, we 
emphasize, as discussed in R99a, R99b, and R99, care must be exercised when interpreting the dis- 
criminative power of these formal limits, since they depend sensitively on the fact that the uniform 
prior has been set to zero outside the range of the parameter space we have considered. 

Figure 2 shows the contours of the two-dimensional posterior distribution for Qrms-PS and Qq, 
derived by marginalizing the four-dimensional distribution over flBh^ and to. These are shown for 
the OVRO and DMR data, for both the open and flat-A models. Constraints on these parameters 
from the OVRO data are consistent with those from the DMR data. 

Figure 3 shows the one-dimensional posterior distribution functions for Clo, CIbH'^, to, and 
Qrms-PS) derived by marginalizing the four-dimensional posterior distribution over the other three 
parameters. From these one-dimensional distributions, the OVRO data favors an open (flat-A) 
model with = 0.33 (0.10), or J^^/i^ = 0.005 (0.005), or to = 10 (11) Gyr, amongst the models 
considered. At 2 cr confidence the OVRO data formally rule out only small regions of parameter 
space. From the one-dimensional distributions of Fig. 3, the data require JIq < 0.69 or > 0.75 (JIq 
< 0.51 or > 0.56), or O^/i^ o.028 (J^b^^ < 0.028), or to < 19 Gyr (to < 19 Gyr) for the open 
(flat-A) model at 2 a. 

While the statistical significance of the constraints on cosmological parameters is not high, it is 
reassuring that the OVRO data favor low-density, young models, consistent with indications from 
most other data. The constraints on Q,Bh? derived from the OVRO data are somewhat puzzling. 
They are more consistent with those derived from the Python and combination CMB anisotropy 
data sets analyzed by R99 and R99b, but less so with those from ARGO (R99a) and more recent 
data sets (Netterfield et al. 2002; Pryke et al. 2002; Stompor et al. 2001) which favor higher 
Ofi/i^. The lower Os/i^ found here is more consistent with the low Cyburt, Fields, & Olive (2001) 
standard nucleosynthesis value determined from helium and lithium abundance measurements, and 
less consistent with the high deuterium-based value of Buries, Nollett, & Turner (2001). 

The peak values of the one-dimensional posterior distributions shown in Fig. 3 are listed in 



See Fig. 4 of R99a, Fig. 2 of R99b, and Fig. 3 of R99 for related cosmological constraints from other data. 
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the figure caption for the case when the four-dimensional posterior distributions are normahzed 
such that -L((5rms-PS = ;uK) = 1. With this normahzation, marginaUzing over the remaining 
parameter the fully marginalized posterior distributions are 1.4 x 10®''(1.3 x 10®'') for the open 
(flat-A) model. This is not inconsistent with the indication from panels a) and b) of Fig. 3 that 
the most-favored open model is marginally more favored than the most-favored flat-A one. 



4. Conclusion 

The OVRO data results derived here are mostly consistent with those derived from the DMR, 
SP94, Python I-III, ARGO, MAX 4-^5, White Dish and SuZIE data. The OVRO data significantly 
constrains Qrms-PS (for the flat bandpower spectrum Qrms-PS = "^^-5 at 1 a) and weakly 
favors low-density, low young models. 

We acknowledge valuable assistance from R. Stompor and helpful discussions with K. Ganga 
and E. Leitch. PM, BR, and TS acknowledge support from NSF CAREER grant AST-9875031. NS 
acknowledges support from the Alexander von Humboldt Foundation and Japanese Grant-in- Aid 
for Science Research Fund No. 14540290. 
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Table 1: Numerical Values for the Zero-Lag Window Function Parameters^ 



360 596 537 753 1.41 



*The value of I where Wi is largest, Zm, the two values of I where Wi = e °'^Wi^, Zg-o.s, the effective multipole, 
k = I{lWi)/I{Wi), and I{Wi) = ESaO + 0.5)Wi/{l{l + 1)}. 



Table 2: Numerical Values for Qrms-PS and 6Ti from Likelihood Analyses Assuming a Flat Band- 
power Spectrum 



Qrms-PS'^ 


Ave. Abs. Err."^ 


Ave. Frac. Err.'^ 




LRd 


(/xK) 






(/xK) 




38 33 


5.5 


14% 


59 

51 


9 X 10^^ 



''The first of the three entries is where the posterior probability density distribution function peaks and the vertical 
pair of numbers are the ±1 a (68.3% highest posterior density) values. 

''Average absolute error on Qrms-ps in /iK. 

''Average fractional error, as a fraction of the central value. 

'^Likelihood ratio. 
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Fig. 1. — Confidence contours and maxima of the OVRO data two-dimensional posterior proba- 
bility density distribution functions, as a function of the two parameters on the axes of each panel 
(derived by marginalizing the four-dimensional posterior distributions over the other two param- 
eters). Dashed lines (crosses) show the contours (maxima) of the open case and solid lines (solid 
circles) show those of the flat-A model. Panel a) shows the (fig/i^, Oq) plane, and panel b) shows 
the {to, 0,0) plane. 
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Fig. 2. — Confidence contours and maxima of the two-dimensional (Qmis-PS) ^o) posterior prob- 
ability density distribution functions. Panel a) shows the flat-A model and panel b) the open 
model. Heavy lines show the ±1 and ±2 a confidence limits and solid circles show the maxima of 
the two-dimensional posterior distributions derived from the OVRO data. Shaded regions show the 
two-dimensional posterior distribution 1 a (denser shading) and 2 a (less dense shading) confidence 
regions for the DMR data (Gorski et al. 1998; Stompor 1997). The DMR results are a composite of 
those from analyses of the two extreme data sets: i) galactic frame with quadrupole included and 
correcting for faint high-latitude galactic emission; and ii) ecliptic frame with quadrupole excluded 
and no other galactic emission correction (Gorski et al. 1998). 



-13- 




Fig. 3. — One-dimensional posterior probability density distribution functions for fio, ^sh^, to, 
and Qrms-PS (derived by marginalizing the four-dimensional one over the other three parameters) 
in the open and flat-A models. These have been renormalized to unity at the peaks. Dotted vertical 
lines show the confidence limits derived from these one-dimensional posterior distributions and solid 
vertical lines in panels g) and h) show the ±1 and ±2 a confidence limits derived by projecting the 
OVRO data four-dimensional posterior distributions. The 2 a DMR (marginalized and projected) 
confidence limits in panels g) and h) are a composite of those from the two extreme DMR data 
sets (see caption of Fig. 2). When the four-dimensional posterior distributions are normalized such 
that ^(Qrms-PS = /xK) = 1, the peak values of the one-dimensional distributions shown in 
panels a) - h) are 1 x 10^^, 2 x 10^'^, 6 x lO^^, 7 x lO^^, 1 x 106^ 2 x 1 x 10^^, and 1 x 106^ 
respectively. 
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